CGL Defect Chaos and Bursts in Hexagonal Rotating non-Boussinesq Convection 



o 
o 



Santiago Madruga" 1 ", Hermann Riecke + , and Werner Pesch* 
+ Engineering Science and Applied Mathematics, Northwestern University, 
Evanston, IL 60208, USA, * Physikalisches Institut, 
Universitat Bayreuth, D-95440 Bayreuth, Germany 
(Dated: February 8, 2008) 

We employ numerical computations of the full Navier-Stokes equations to investigate non- 
Boussinesq convection in a rotating system using water as a working fluid. We identify two regimes. 
For weak non-Boussinesq effects the Hopf bifurcation from steady to oscillating (whirling) hexagons 
is supercritical and typical states exhibit defect chaos that is systematically described by the cubic 
complex Ginzburg-Landau equation. For stronger non-Boussinesq effects the Hopf bifurcation be- 
comes subcritical and the oscillations exhibit localized chaotic bursting, which can be modeled by a 
quintic complex Ginzburg-Landau equation. 

PACS numbers: 47.52.+j,5.45. Jn,47.20.Lz,47.27.Te 
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The complex Ginzburg-Landau equation (CGL) as the 
universal description of weakly nonlinear oscillations has 
been studied theoretically in great detail. The classical, 
supercritical case, which involves only two independent 
parameters, looks deceptively simple. It exhibits, how- 
ever, a vast variety of qualitatively different phenomena 
including different types of spatio-temporally chaos (e.g. 
[j], 0). If the oscillations arise in a subcritical Hopf bi- 
furcation the quintic CGL comes into play, which intro- 
duces further interesting states including, for instance, 
intermittent bursts in the oscillation amplitude. 

In contrast to the extensive theoretical work on 
the CGL, direct experimental validation of its various 
regimes of complex behavior are rather scarce, in partic- 
ular for the two-dimensional case 3]. In this paper we 
present detailed numerical computations for an experi- 
mentally realizable thermal convection experiment that 
exhibit defect chaos and bursts. We show that for the 
defect chaos the cubic CGL should provide a systematic 
description, while the bursts can be modeled with a quin- 
tic CGL. 

Rayleigh-Benard convection of a fluid layer heated 
from below in systems with large aspect ratio, in which 
the lateral dimension L of the layer is much larger than 
its thickness h, has proved to be a paradigmatic experi- 
mental system 4] for studies of complex patterns. Above 
a critical temperature difference AT C across the layer, 
which corresponds to the critical value R c of the dimen- 
sionless Rayleigh number one observes in the simplest 
case the familiar striped (roll) patterns with wavenum- 
bers q close to the critical wavenumber q c . However, in 
systems in which AT C is large fluid properties like the 
thermal expansion coefficient or the viscosity vary sig- 
nificantly across the layer. Under these non-Boussinesq 
(NB) conditions the instability of the homogeneous state 
leads to hexagonal convection patterns [5| . 

If the chiral symmetry of the system is broken by rotat- 
ing the layer about a vertical axis with frequency new 
interesting dynamics arise. For Q above the Kiippers- 
Lortz frequency VLrl one finds in the Boussinesq case 
immediately at onset R c domain chaos in which rolls per- 
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FIG. 1: Sketch of bifurcation diagram for rotating NB- 
convection. Solid (dashed) lines denote stable (unstable) 
branches. For the oscillating ('whirling') hexagons maxima 
and minima of the oscillations are indicated (thick solid lines) . 
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FIG. 2: Disordered state of whirling hexagons in case A for 
e = 0.2 (R = 7396.6), L = 32 • 2n/q c = 38.3. a) Snapshot of 
0(r). b) Magnitude \7i\ of oscillation amplitude of snapshot 
shown in a). Red (blue) indicates large (small) values of \7i\. 



sistently switch orientation (cf. 3). In contrast, in the 
NB case the hexagons are steady in this regime. How- 
ever, weakly nonlinear theory predicts a secondary oscil- 
latory instability at e = (R — R c )/R c = e# that leads 
to 'whirling hexagons' in which the three hexagon am- 
plitudes oscillate about their mean, with their phases 
shifted by 2tt/3 relative to each other @, 0, H, (cf. 
FigQ). 

Here we perform computer experiments on rotating 
non-Boussinesq convection with water as working fluid. 
The experimental realization of these simulations should 
pose no serious difficulty in presently available experi- 
mental setups. Focusing on the whirling hexagons, we 
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obtain different spatio-temporally disordered states, an 
example of which is shown in Fig|2t> and use the complex 
oscillation amplitude TL (see (0 below) to interpret them 
in terms of suitable complex Ginzburg-Landau equations. 
The magnitude \TL\ is shown in FigJSjD. 

Specifically, we consider two situations correspond- 
ing to weak and intermediate NB effects, respectively, 
and solve the full Navier-Stokes (NS) equations with 
leadin g-or der NB effects included, as discussed previ- 
ously [TTj | . We obtain spatially periodic stationary 
hexagon patterns and test their stability using a standard 
Galerkin method. The temporal evolution of the states 
is then simulated with a pseudo-spectral code with pe- 
riodic lateral boundary conditions [HI E[ • We consider 
two set-ups differing in the mean temperature To of the 
fluid layer. In case A we choose To = 14° C resulting in 
AT C = 6.4°C and a value of Q = -2.2 for Busse's NB 
parameter 5]; case B with T = 12°C and AT C = 8.3°C 
is more non-Boussinesq and yields Q = —3.62 [2^. In 
both cases h = 0.492 cm and Q = 65 v$/h 2 with vq being 
the viscosity of the fluid in the mid-plane. 

In case A we obtain above the Hopf bifurcation at e# 
states in which almost all convection cells are oscillating 
about a mean, but for general initial conditions the os- 
cillations of different cells are out of phase with respect 
to each other. This is shown in Fig|2ji where the solution 
is visualized by the deviation #(r, t) of the temperature 
from the conduction profile in the mid-plane of the fluid 
layer. Here r denotes the horizontal coordinates. To ex- 
tract the oscillations explicitly we make use of the fact 
that the underlying hexagon pattern itself is well ordered. 
We first demodulate each snapshot in space by writing 



<9(r, t)=J2 A o( r i l ) eiqcilyV + c - c - + h -°- t - 



(1) 



Here the wavevectors q c hj represent the three dominant 
wavevectors of the hexagons and h.o.t. denotes their har- 
monics. For strictly periodic stationary hexagon patterns 
the amplitudes Aj, j = 1..3, do not depend on r or t and 
have the same modulus Ah ex - They can be chosen real 
and positive in our case (water). More generally, the am- 
plitudes Aj (r, t) depend slowly on space representing the 
contributions from the side-bands of the basic wavevec- 
tors q c hj. 

Slightly above the Hopf bifurcation the amplitudes Aj 
vary in time. For patterns that are close to periodic in 
space they can be expressed as 

Aj = (A hex + /i 1 / 2 [e^^ 3 He iujt + c.c] + 0(/x)) x 

exp (iAqhj • r + i/i 1/2 Mr,t)), (2) 

where \i = (R — Rh)/Rh and uo is the Hopf frequency. 
Here TL denotes the complex oscillation amplitude, which 
we extract by standard demodulation of a series of snap- 
shots like Fig|2K in space and time. The space depen- 
dence of the phases <j)j captures slight deformations of 
the underlying hexagon lattice. 
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FIG. 3: a) Bifurcation diagram for the subcritical case B 
(cf. FigQwith the rolls omitted), b) Bursting state obtained 
with NS-simulation in case B for e = 0.5 (R = 6121.9) and 
L = 16 • 27r/q c = 19.1. \TL\ color-coded as in FigEJ), contour 
lines indicate hexagon pattern. 



A representative snapshot of the magnitude of the os- 
cillation amplitude is shown in FigJ2t>. A more detailed 
analysis shows that the domains with low oscillation mag- 
nitude are due to spiral defects in the complex oscillation 
amplitude, i.e. TL = at isolated points. During the evo- 
lution of the system these defects are persistently created 
and annihilated in pairs. In appearance, this state is thus 
quite similar to the defect chaos of the cubic CGL (see 
below) . 

When the NB effects are stronger ( case B ) the Hopf 
bifurcation is shifted to larger values of e and, unexpect- 
edly, becomes subcritical. The numerically determined 
bifurcation diagram showing the jump in the oscillation 
amplitude, the small hysteresis, and the restabilization 
(reentrance [ll|) of the steady hexagons is presented in 
FigEJi. As a consequence, the dynamical states obtained 
in this regime are very different, as shown in FiglHt). 
The contour lines indicate the underlying hexagon pat- 
tern and, as in Fig|2t), the colors represent the oscillation 
magnitude. Strikingly, the oscillations are now localized 
into relatively small domains. Only within these domains 
the convection cells are elongated to footballs and are 
whirling. In most parts of the system the hexagons are 
relatively steady. The burst-like temporal evolution of 
the hexagon amplitudes is evident from Fig^I It shows 
the normalized oscillation intensity in the bursts, which 
we define as I(t) = TV -1 J^ >Q 5 \ n \ \Tl\ dxdy where TV 
is the temporal mean of the integral. The intensity I(t) 
exhibits substantial intermittency reflecting the growth 
and decay of bursts (thick solid line), while in case A 
(thin solid line) the intensity fluctuates relatively little. 

To provide insight into the complex dynamics obtained 
in our numerical experiments we make use of the fact that 
the whirling hexagons arise from a Hopf bifurcation. It 
is therefore natural to expect that in the supercritical 
case the dynamics can be captured by the CGL for the 
oscillation amplitude, 



d T H = fJL 2 n + dV 2 H - cH\H\* 



(3) 



We determine the linear coefficients \±2 and d = d r + idi 
using the Galerkin stability computations and extract the 
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FIG. 4: Normalized oscillation intensity I(t) (see text). Case 
A (thin black line) and case B (thick red line). 




FIG. 5: Bursting state obtained with quintic CGL for param- 
eters corresponding to FigEI M2 = 0.064, d = 1.90 + 0.033i, 
c — — 1.10 + 7.19i, g — 3.61 + 1.46z, except for the larger 
system size, L — 100. a) magnitude \Tt\ of oscillation, b) en- 
larged portion from a): magnitude \TL\ (top) and associated 
wavenumber |k| of oscillation (bottom). 



nonlinear coefficient c = c r + ici from the dependence 
of the oscillation frequency on the oscillation amplitude. 
We find d = 1.3 + 0.05H and c = 0.95 + 9. 16i. Consistent 
with predictions based on a theory for small-amplitude 
hexagons 0, |?j , these values position the system in the 
Benjamin- Feir- stable regime (l + diCi/d r c r > 0) in which, 
however, defect chaos persists due to the wavenumber 
selection by the spiral defects yj,|2||. Indeed, we find in the 
NS-simulations that spatially homogeneous oscillations 
are stable. To reach the chaotic attractor shown in Fig0 
we employ initial conditions in which the oscillations are 
not synchronized across the system. 

To capture the bursting dynamics found in the subcrit- 
ical case B the CGL ((3} has to be extended to include 
a quintic term, — g\H\ A TL. Proceeding as in case A, we 
determine the coefficients of the resulting quintic CGL 
and simulate it in two dimensions. As shown in FigEt, 
which depicts a snapshot of the magnitude |W|, we obtain 
quite similar bursts. Based on this qualitative agreement 
we use the CGL to gain further insight into the bursting 
found in the NS-simulations. 

The burst mechanism can be elucidated by consider- 
ing the local gradient of the phase of H = IZe 1 ^ , i.e. the 
wavevector k = of the oscillations. Note that k is 
not related to the wavevectors making up the hexagon 
pattern. Gradients in \TL\ induce significant differential 
phase winding due to the strong amplitude dependence 
of the oscillation frequency (\ci/c r \ ^> 1). This leads to a 
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FIG. 6: a) V • <fi generated by burst shown in Fig|3j3. b) 
Growth rate a of the oscillatory mode as a function of the 
wavenumber q = q c + Aq. 



buildup of |k| at the perimeter of the bursts while |k| re- 
mains small at their cores, as illustrated in FigEb, which 
shows an enlargement of a burst (top) and the associated 
|k| (bottom). The increased wavenumber k enhances the 
dissipation via the diffusion term in © and eventually 
leads to a collapse of the oscillation amplitude. This 
mechanism has been studied previously in one dimension 
|l2l Il3l \l4 . Il5l llvL Il7j . It has been shown to underly 
the dispersive chaos observed in binary-mixture convec- 
tion [13, 13 anci , interestingly, can be strong enough to 
avoid blow-up even if there is no saturating nonlinear 
term at all (c r = g r = 0) 0,0,0. The collapse can 
also be interpreted in terms of two colliding fronts that 
connect the steady base state with the nonlinear oscil- 
latory state. Since these fronts select a non- vanishing 
wavenumber for the nonlinear state, the base state can 
invade the nonlinear state even above the Hopf bifurca- 
tion ('retracting fronts') 0. Note that the collapse of 
the bursts described here is not due to a break-down of 
the underlying hexagonal structure 0, . 

Assessing the validity of the CGL (0 and its quintic 
extension it is clear that for such a secondary Hopf bifur- 
cation a systematic treatment must also include the cou- 
pling of the oscillation mode to slow, long-wave deforma- 
tions of the hexagon lattice [2^. They are characterized 
by gradients in the phase vector <f> = (d>i, (02 — 4>s)/V3j 
|8|- In particular, local deformations of the hexagon lat- 
tice enter through an additional term ~ which 
captures the dependence of the growth rate a of the oscil- 
lations on local compressions of the lattice. Conversely, 
gradients in Ti drive <fi |8| . As a careful analysis of the NS- 
simulations shows, they lead to a decrease of V • <fi inside 
the burst, as shown in FigEJi. This has the same effect 
on a as a decrease in the wavenumber of the hexagons, 
Qc —> Qc + with Aq — V • 0/2, and therefore results in 
an increase in a (see FiglS^), providing a positive feed- 
back that enhances the bursting activity compared to the 
quintic CGL. 

Of course, the coupling of the oscillations to the under- 
lying lattice exists also in the supercritical case shown in 
FigEl Consistent with the weakly nonlinear description 
of rotating hexagons H,0, we find, however, that in this 
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regime the growth rate a depends only weakly on the lo- 
cal wavenumber of the hexagons (see FigEb) and varies 
by less than 10% across the system. We therefore sug- 
gest that the classic cubic CGL should describe whirling 
hexagons very well in the supercritical case A. 

In conclusion, we have numerically investigated 
whirling hexagons that arise in rotating NB-convection 
and have identified two different regimes. In the weakly 
NB case the oscillations typically exhibit defect chaos and 
our analysis suggests that this state should be well de- 
scribed by the cubic CGL. Rotating NB-convection would 
then represent one of only a few experimentally realizable 
physical systems |2| in which at least one of the complex 
states of the two-dimensional CGL is accessible. 

For stronger NB effects we found that the Hopf bi- 
furcation becomes subcritical and typical states exhibit 
bursts. Such bursts and the related retracting fronts have 
been discussed in some detail previously in one dimen- 
sion [H [H El [H EH El. The intermittent bursting 
behavior, however, has not been studied in detail yet. In 
two dimensions even less is known. For instance, the con- 
ditions for the persistence of bursting when stable steady 
and whirling hexagons compete remain to be understood. 
Preliminary simulations indicate that for e = 0.7 the 
bursting persists but not for e = 0.9 (cf. FigJSK)- 

Since the oscillations arise in a secondary bifurcation 
the oscillatory mode is, in principle, coupled to the defor- 
mations of the hexagon lattice. While in the supercritical 
case our results indicate that this effect should be small, 
it is more significant in the subcritical case, where de- 
formations of the lattice modify the growth rates of the 
oscillatory mode substantially (cf. FigEb). Compared 
to the bursting behavior of the quint ic CGL alone, this 
leads to an enrichment of the scenario, which requires 



further investigations. 

In our simulations we have taken great care to obtain 
defect-free hexagonal lattices. Since boundaries tend to 
introduce defects we have performed some runs in which 
circular ramps in the Rayleigh number mimic a circu- 
lar container. By stabilizing the hexagon pattern near 
the boundary through an additional patterned volume 
heating we were able to recover the defect chaos and the 
bursting state. Conversely, it would be interesting to al- 
low the underlying lattice to be disordered. In that case 
new phenomena may arise from the interaction of the 
penta-hepta defects with the whirling mode and signifi- 
cantly more complex states may be found (cf. 10]). 

While our numerical experiments have been very pro- 
ductive in identifying and exploring various complex 
states in rotating NB-convection, current computational 
limitations do not allow us to use the NS-simulations 
to investigate in detail the statistical properties of these 
states. For instance, for the defect chaos we cannot 
address the possibility of deviations from the squared 
Poisson distribution for the defect statistics (2£| or the 
expected transition to glassy states with exceedingly 
slow dynamics |2l| as the wavenumber of the underlying 
hexagon pattern is changed Similarly, investigations 
of the ramifications of the hexagonal anisotropy due to 
the underlying lattice or of the broken chiral symmetry 
may require system sizes that are still beyond current 
computational capabilities. Questions like these can so 
far only be investigated in experiments. 

Two of us (HR, WP) dedicate this paper to the mem- 
ory of our good friend Lorenz Kramer. 

We gratefully acknowledge support from the Depart- 
ment of Energy (DE-FG02-92ER14303). 
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